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Even though the Dreh-Yan-West formulation is the most rigorous and 
weh-established framework to compute the exclusive processes, its utility has 



Q^ \ been limited only to the spacelike region because of the intrinsic kinematic 



\^ \ constraint g+ = 0. We present an explicit example demonstrating how one 

fvq ! may obtain the necessary information (i.e. nonvalence or so called Z-graph 

(N , 

^O i contribution) in the timelike region of exclusive process without encountering 



a formidable task of direct calculation that has hindered so far the progress 
in this area. In the analysis of qQ bound state form factors using an exactly 



Qh! solvable model of (3+1) dimensional scalar field theory interacting with gauge 

fields, the results analytically continued from the spacelike region coincide 



k> \ exactly with the direct results in the timelike region. This example verifies 

H ; 

5t , that the method of analytic continuation is capable of yielding the effect 

of complicate nonvalence contributions. The meson peaks analogous to the 
vector meson dominance(VMD) phenomena are also generated at the usual 
VMD positions. 
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I. INTRODUCTION 

The Drell-Yan-West (g+ = g° + g^ = 0) frame in the hght-front quantization provided an 
effective formulation for the calculation of various form factors in the spacelike momentum 
transfer region q^ = —Q^ < [l|]. As an example, only parton-number-conserving Fock 
state (valence) contribution (e.g. Fig. 1(a)) is needed in g"*" = frame when the "good" 
components of the current, j+ and j_l = {jx,jy), are used [0 for the spacelike electromagnetic 
form factor calculation of pseudoscalar mesons. Successful light-front quark model (LFQM) 
description of various pseudoscalar form factors can also be found in the hteratures |^|-[7|. 
However, not all is well. For the higher spin system composed of fermions, one needs 
in general to consider the nonvalence contribution (e.g. Fig. 1(b)) to fulfill the covariance 
requirement even when j"*" and j± are used in the g+ = frame. The nonvalence contribution 
in the g+ = frame is called the "zero-mode" . The details of zero-mode contributions for 
the spin-one system with the j~^ current coupled with longitudinal polarization vector as 
well as the semileptonic transition form factor of pseudoscalar meson with the j± current 
can be seen in Refs. [||-01- 

Furthermore, the timelike (g^ > 0) form factor analysis in the light-front quark model 
has been hindered by the fact that g"*" = frame is defined only in the spacelike region 
(g^ = g^g^ ~ Qi. < 0)- While the g^ 7^ frame can be used in principle to compute the 
timelike form factors, it is inevitable (if g^ 7^ 0) to encounter the nonvalence diagram arising 
from the quark-antiquark pair creation (so called "Z-graph") as shown in Fig. 1(b). The 
main source of difficulty in calculating the nonvalence diagram (Fig. 1(b)) is the lack of 
information on the black blob which should contrast with the white blob representing the 
usual light- front valence wave function. In fact, it was reported that the omission of 
nonvalence contribution leads to a large deviation from the full results. The timelike form 
factors associated with the hadron pair productions in e~^e~ annihilations also involve the 
nonvalence contributions. Therefore, it is not a simple task to explore the timelike region 
with the framework at hand and the investigation is still just at the beginning stage to the 



best of our knowledge. 
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FIG. 1. The light-front quark model description of a electroweak meson form factor: (a) the 
usual light-front valence diagram and (b) the nonvalence(pair-creation) diagram. The vertical 
dashed line in (b) indicates the energy-denominator for the nonvalence contributions. While the 
white blob represents the usual light-front valence wave function, the modeling of black blob has 
not yet been made. 

The aim of this paper is to provide at least a clear example of demonstration showing 
how one might take advantage of the existing formulation to get the new information (i.e. 
nonvalence contribution or so called Z-graph) necessary in the timelike region of exclusive 
processes. This can be done by the analytic continuation from the spacelike form factor 
calculated in the Drell-Yan-West (g^ = 0) frame to the timelike region. In this paper, 
we focus on the case that does not require the zero-mode contributions even though one 
can in principle apply the same method presented in this work to the case of zero-mode 
contributions once they are identified. However, in the case of the zero-mode contributions, 
an additional step is necessary to identify themQ. 



^An example including the zero-mode contribution can be shown in the calculation of the /_ form 
factor for the exclusive semileptonic decays between two pseudoscalar mesons using the fermionic 
covariant model with the monopole smearing vertices [@,^. In this case, we indeed verified the 
applicability of our method by checking if the identical results are obtained when different com- 



For an explicit demonstration of our work, we use an exactly solvable model of (3 + 1) 
dimensional scalar field theory interacting with gauge fields. The Lagrangian of the system 
is given by ||rT| 



+ (a^0fe + iebA^(j)by{d^(f)b + iebA^(j)b) - ml4(f)b 

_1(M2$2 _ a^$a^$) + ^$(0t0, + 0t^J, (1) 

where 4>a{b) corresponds to the bosonic quark field with the mass ma(6) and $ is the bound 
state meson with the mass M. Our model is essentially the (3 + 1) dimensional extension 



of Mankiewicz and Sawicki's (1 + 1) dimensional quantum field theory model ||T^, which 
was later reinvestigated by several others pn] , p!3| -[T6[] . The starting model wave function is 
the solution of covariant Bethe-Salpeter (BS) equation in the ladder approximation with a 
relativistic version of the contact interaction |]T^. Here, we do not take the Hamiltonian 
approach. The covariant model wave function is a product of two free single particle prop- 
agators, the overall momentum-conservation Dirac delta, and a constant vertex function. 
Consequently, all our form factor calculations are nothing but various ways of evaluating 
the Feynman perturbation-theory triangle diagram in scalar field theory. 

Even though the model may not be as realistic as one may want, it achieves the goal 
of being the first example up to our knowledge. We think that it provides at least some 
guidance how one can persue the same (or similar) idea in the more realistic (perhaps more 
phenomenological in that respect) models. In this model, we calculate: (A) the timelike 
process oi j* —>■ M + M transition in g+ 7^ (g^ > 0) frame, (B) the spacelike process of 
M ^ 7* + M in g+ 7^ (g^ < 0) frame, and (C) the spacelike process of M ^ 7* + M in 
g+ = frame. Using the analytic continuation from g^ < to g^ > 0, we show that the 
result in (C), which is obtained without encountering the nonvalence contributions at all. 



ponents of the current are used in different frames, i.e., q~^ j^ and g^ = frames. See also the 
discussion in Section V. 



exactly reproduces the result in (A). In fact, all three results (A), (B), and (C) coincide with 
each other in the entire q^ range. We also confirm that our results are consistent with the 



dispersion relations |T7|-pO[]. We consider not only for the equal quark/antiquark mass case 
such as the pion but also for the unequal mass cases such as K and D. 

The paper is organized as follows: In Sec. II, we derive the timelike electromagnetic (EM) 
form factor of 7* — ;> M + M process in the g^ 7^ frame (A) and discuss the singularities 
occuring from the on-energy shell of quark-antiquark pair creation. In Sec. Ill, the spacelike 
form factor of M — i> 7* + M process is calculated both in the q^ ^ Q (B) and g^ = (C) 
frames. We then analytically continue the spacelike form factors to the timelike region. The 
singularities occured in the timelike region are also discussed. In Sec. IV, for the numerical 
calculation of the EM (vr, K, and D) meson form factors for three different cases (A), (B) and 
(C), we use the constituent quark and antiquark masses {niu = Tn^ = 0.25 GeV, mg = 0.48 
GeV, and rric = 1.8 GeV) []7|,|3|,|21|| and show that the form factors obtained from those three 
different cases are indeed equal to each other for the entire g^ region. The meson peaks 
analogous to the vector meson dominance (VMD) are also obtained. The conclusion and 
discussion follows in Sec. V. 

II. FORM FACTORS IN THE TIMELIKE REGION 

The EM local current J^(0) responsible for a virtual photon decay into two qQ bound 
states in the scalar field theory can be calculated using the diagrams shown in Fig. 2. The 
covariant diagram shown in Fig. 2(a) is equivalent to the sum of two light-front time-ordered 
diagrams in Figs. 2(b) and 2(c). The EM current J^(0) obtained from the covariant diagram 
of Fig. 2(a) is given by 

J^(0) = ze,/ fd'k- , . / „ . (g - 2kr ^ 

J [q — ky — m^ + ze 

1 



(g — kY — ml + ie {q ~ k — P2Y — m^Q + ^^ 



^7^ 2 ~ ~'~ ^^(mg <-»• mg of the first term), (2) 

where Trigi^Q^ and eg(Q) are the constituent quark (antiquark) mass and charge, respectively. 
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FIG. 2. The electromagnetic decays of a photon into two-body bound sates, i.e., 7* — s- qq{ox 
QQ)—f M.{qQ)M.{qQ), in scalar theory: Covariant representation (a), and the light-front time 
ordered contributions to the decay amplitude (b) and (c). 

The corresponding form factor F{q^) of the qQ bound state in tinielike(g^ > 0) region is 
defined by 



J^(0) = (Pi - P2rF{q' 



(3) 



where q = Pi + P2, Pf = P2 = M"^ and M is the mass of a qQ bound state scalar particle. 
Using the Cauchy integration over k~ in Eq. (0), we can find each time-ordered con- 
tribution (Figs. 2(b) and 2(c)) to the timelike form factor -F(g^ > 4M^) in Eq. (^. This 
procedure allows us to analyze the singularity structure of each light-front time-ordered di- 
agram as well. For the calculation of each light-front time-ordered contribution, we take the 
purely longitudinal momentum frame, i.e., q^ ^ 0, q± = and Pi± = P2± = 0. Accordingly, 
g^(= q^q~) > ^M"^ is given by 



g2 = M2(l + a)Va, 



(4) 



where a = P2 / Pi = q^ / Pi —1 is the longitudinal momentum fraction and the two solutions 
for a are given by 



a± 



(J^ 



1 ± 



7 g2 



1-1. 



(5) 



V2M2 ; V V2M2 

Note that both a± = 1 correspond to the threshold, q^ = 4M^. The EM form factor -F(g^ 
in Eq. (^ is independent of the subscript sign of a. Thus, one can take either a+ or a. 
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to calculate F{q^). Here, for convenience, we use a = a^ which ranges from to 1 for the 
physical momentum transfer region, i.e., a_ ^ as g^ ^ oo and 1 as g^ -^ 4M^. Of course, 
we can use a = a+ equally well and verify that the two results (a+ and a_) are exactly 
same for the calculation of F{q^). 

Since g+ > Pi > P2 > for a = a_, the Cauchy integration over k~ in Eq. (^ 
has two nonzero contributions to the residue calculations, one coming from the interval (i) 
< k+ < P^ (see Fig. 2(b)) and the other from (ii) P^ < k^ < q^ (see Fig. 2(c)). 
The internal momentum k~^ is defined by k'^ = xP^, where x is the Lorentz invariant 
longitudinal momentum variable. The "good" component of the current, J~^{0), is used in 
our computation of the two light-front diagrams Figs. 2(b) and 2(c). In the following, for 
simplicity, we will explicitly write neither the obvious second term in Eq. (0) nor the charge 
factor (cg or eg). 

In the region oi < k'^ < Pi~, the residue is at the pole of k~ = [m^ + k\ — ie]/k~^, 
which is placed in the lower half of complex-Zc" plane. Thus, the Cauchy integration of J^ 
in Eq. (]^) over k~ in this region gives 

'■^i^ ,.^ ,0, Q+ -2k- 






1 



{ml + kl)lk+ - {ml + kl)l{q+ - k+) 



X-F 



T, (6) 

q- - P^ - {ml + kl)/k+ - (m| + kl)l{P+ - k+) 

where the subscript h in J^(0) implies the current in Fig. 2(b) and g is fixed by the 
normalization of the EM form factor at g^=0. The last two terms in Eq. (^ correspond to 
the two and three particle intermediate states. We represent the energy denominators of the 
two and three particle intermediate states as E\ and E\ in Fig. 2(b), respectively. 
To analyze the singularities of Eq. (^), we further integrate over k^ and obtain 
,, , , 99/"^, x{l + a-2x)/{l + a) ^ f £^\ 

j*(o) = .y/^d.J — J — ii„(^), (7) 

where £l = x{l + a — x)M'^/a — m? and £l = x{l — x)M^ — [xm^ + (1 — 2;)m^]. While S^ is 



not zero (£^| 7^ 0) in general for the entire physical region, S2 can be zero when g^ > Am?,Qs. 
The singular structure of £\ — £\ term in Eq. (0) depends on whether a qQ bound state 
scalar particle is strongly bounded (M^ < m^ + toq) or weakly bounded (M^ > ut^^ + toq). 
As we will show in our numerical calculations (Section IV), anomalous threshold appears for 
M^ > m? + m?0 while only the normal threshold of bound state exists for M^ < m? + ttlq. 



In the region of Pi < k^ < q^, the residue is at the pole oi k = q — [m?„ + {q^ — k 



-l; 



ie]/(g^ — A;+), which is placed in the upper half of complex-fc plane. Thus, the Cauchy 
integration of ./"^(O) in Eq. (H) over k~ in this region yields the result 

2k+ 



J+(0) = -Tig' / dk+d'k_ 



q 

k+{q+ -k+){P^ -k+) 
1 



q' 



K + kl)lk+ - K + A;i)/(g+ - A;+) 



1 



T' (S) 

g- - P{ + (m| + kl)/{Pt - k+) - {ml + kl)/{q+ - k+) 

where the subscript c in J^i^) means the current in Fig. 2(c). After the integration over 
the k± in Eq. (Rl), we obtain 



J^o) 



2.2 f\^aX{l + a-2aX)/{l + a)_^jS, 



-TT'g' I dX 

'0 



Cc Cc 



£!i 



(9) 



where x = 1 + a{l - X), S^ = X(l - X)M'^ - [Xm| + (1 - X)ml] and £| = aX{l + a- 
aX)M'^/a — m^. The pole structure in Eq. (^ is equivalent to that of Eq. (^. 
Consequently, the timelike form factor in Eq. (|^) is given by 



F{q' 



n'g' 



a^ 



IJo 



1 rx(l + a-2x), /£. 
dx{ -7 — ^In' ■ 



cb 



£'0 



cb 



ax{l + a — 2ax) / £2 

cc cc \ cc 

^3 ~ ^2 \^3 



(10) 



where £^21^3)^2 ^'^d £^ are defined in Eqs. (|^ and (y). 



III. FORM FACTORS IN SPACELIKE REGION AND THE ANALYTIC 
CONTINUATION TO THE TIMELIKE REGION 



In this section, we calculate the EM form factor in spacelike momentum transfer region 
and then analytically continue to the timelike region to compare the result with the timelike 
form factor (i.e. Eq. ([T0|)) that we obtained in the previous section. The EM current of a 
qQ bound state in spacelike momentum transfer region is defined by the local current j^(0); 



j'^(0) = (Pi + P2)^^(g' 



(11) 



where q = Pi — P2, q^ < and T{q^) is the spacelike form factor. The EM current j'^(O) 
obtained from the covariant triangle diagram of Fig. 3(a) is given by 

1 _ „ _..„ 1 1 



j^(0) = ig' d'^k 



{Pi — ky — ml + ie 



(Pi + P2- 2ky 



(P2 ~ ky — m? + iek'^ — uiq + ie 



(12) 



Q 






(b) 



(c) 



FIG. 3. Covariant triangle diagram (a) is represented as the sum of light-front triangle diagram 
(b) and the light-front pair-creation diagram (c). 

As in the case of the timelike form factor in Sec. II, the Cauchy integration of k~ in 
Eq. ([T^) has also two contributions to the residue calculations, one coming from the interval 
< A;^ < P2 (see Fig. 3(b)) and the other from P2 < k^ < Pi (see Fig. 3(c)). Again, only 
the "good" -current j^(0) in Eq. ( 1T2|) is used to obtain the contributions from Figs. 3(b) and 
3(c). 



In the region of < A;+ < P2 , the residue is at the pole of k~ = [171% + k'j_ — ie]/k~^, 
which is placed in the lower half of coniplex-/c~ plane. Thus, the Cauchy integration of j~^ 
in Eq. (jT^) over k~ in this region yields 

J, (U) - ng J^ dk d k^ ^^^ _ ^^^^^^ _ ^ _ ^^^^^ 



X-F 




(13) 
P^- -q-- [ml + {k^ + q±y]/{P^ - k+) - [m| + kl]/k+ 

where the subscript h in j^(0) implies the current in Fig. 3(b). 

In the region of P2 < k^ < P^, the residue is at the pole of k^ = P{ — [im? + k\ — 
it\/{Pi — k^), which is placed in the upper half of complex-/;;" plane. Thus, the Cauchy 
integration of j^ in Eq. ([12D over k~ in this region becomes 

(Pi + P2- 2k)+ 



Jc (0) = ng' 



p+ 



dk^d^k I 



(Pi - A;)+(Pi -k- q)+k^ 
1 



P{ - [ml + fci]/(P+ - k+) - [ml 



kl]/k^ 



P2 


- P,- + K + kl]/iP^'- - 


-k+)- 


-K + (g^ + A;^)2]/(P+- 


-k+) 



(14) 



As one can see from Eq. (0), the nonvalence contribution (Fig. 3(c)) vanishes only in 
the g''" = frame. Note that, for the spin-one system composed of fermions, however. 
Equation (jT^) does not vanish even in g+ = frame when the longitudinal polarization 
vector is involved as mentioned in the introduction. In the following, we investigate the 
spacelike form factor J^{q^) given in Eq. (0) using both g^ 7^ and g^ = frames. We 
then analytically continue to the timelike region in order to compare the result with the 
direct calculation of the timelike form factor F{q'^) presented in the previous section. 
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A. The purely longitudinal(g"'" 7^ and g^=0) frame 

In the purely longitudinal momentum frame q^ 7^ 0, g^ = 0, and Pix = P2L = 0, 
the momentum transfer q^ = q^q~ can be written in terms of the longitudinal momentum 
fraction (3 = P^/P^ = 1 - q^/P^] 

q^ = -M^{l-l3f/l3<Q, (15) 

where the two solutions of (3 are given by 



P±={^-47,)^\{^-4h\-^- (16) 



2M2y V V 2M2y 

The form factor J-'{q'^) in Eq. (0) is also independent of the subscript sign of /5. However, the 
condition P^ < P^ was used in obtaining Eqs. (|13]) and ([T^) and thus here we use (3 = (3- 
(0 < /3 < 1) in spacelike region. As shown in Eqs. (|T2|)-(|l^, the sum of valence (Fig. 3(b)) 
and nonvalence (Fig. 3(c)) diagrams is equivalent to the covariant triangle diagram in Fig. 
3(a). 

For the analysis of singularity structures, we integrate over k± and obtain from the 
valence contribution(Fig. 3(b)); 



'0 £|-£2 Vfj 



g{0) = rcV d. ^-^'^'Z^.r' mr^]. (17) 



where £^ = x{l — x)M'^ — [xm'^ + (1 — a;)m^] and ^| = Px{l — Px)M'^ — [Pxm'^ + (1 — Px)m%]. 
It turns out that Eq. (P!7| ) has no singularities because £^| 7^ 0, £^2 7^ O5 ^^d £^| — £^2 7^ ^^ 
the entire g^ region. On the other hand, the kj_ integration for the current j^{0) in Eq. (H) 
yields 

,^(0)^.V /'.;>- '^-^;^'^^-^ n(gl (18) 



where a; = 1-(1-/5)A', S^ = {l-(3)X[l-{l-[3)X]M^-[{l-f3)Xml + {l-{l-(3)X)ml] and 
£^| = —{1 — j3)'^X{l — X)M'^ / (3 — m? . While £| corresponding to the energy denominator of 
the two particle intermediate state does not vanish both in spacelike and timelike region, £^ 
of the three-particle intermediate state is not zero only for the spacelike momentum transfer 
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region. For the timelike region, £^ can be zero so that the singularities start at g^j^ = Am^Q. 
for 'y*qq{'y*QQ) vertex. The singularity structure of £^| — S^ in Eq. (|18]) is the same as in 
the case of timelike form factor (Section II), following the condition of a qQ bound state. 

The EM form factor J^{q^) in Eq. (^) of a qQ bound state in g+ 7^ frame is then 
obtained by 



I + 13 Jo I S^-8l \£l 



^(g^g+^0) = -^ dx\ " ^ - " - " M n 



l-(3yx{2x-l)^ (8, 



where 82,8!^, £2 ^^^ ^2 ^^^ defined in Eqs. (|1^) and ([T8|). Here, /? is a function of g^. 
According to the analytic continuation, the sign of g^ in Eq. ([T9|) must be changed from — 
to + for the timelike region. 

B. The Drell-Yan-West (g+=0) frame 

In Drell-Yan-West frame, q~^ = 0, q^ = —q'j_, and Pi± = P2± = 0, the '+' component of 
the current has only the valence contribution, i.e., j^{0) in Eq. (|T3|). The current jfl'iO) in 
g+ = frame is given by 

where i± = k^ + xq±/2, A = x(l — x)M^ — [xm^ + (1 — x)m|], and ^^ = x^g^/4. The angle 
(0 < < 27r) is defined by i'j. ■ q_L = |^±||q±| cos0. 
Integrating Eq. (|20D over and i±, we obtain 



,.(o).-^^i!£L^/',. (1- 



y'gi(4M2 + g^)^o y^(x-x+)(x-x_) 
X tanh~^ 



q\ X 



[\l 4M2 + ql ^[x-x+){x-x^) 



(21) 



where 



2(M2 -ml + m%) 



4(M2 - m2 -K m^)2 Avfij 



4M2 + g2_ ^ (4M2 + g2)2 (4M2 + gi)' 
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(22) 



The analytic continuation from spacelike to timelike region in the g"*" = frame requires the 
change of q± to ig_L in Eqs. (|21| ) and (^). We note from Eqs. ( pT]) and ( p2D that the result of 
the timelike region exhibits the same singularity structure as the direct analyses in g+ 7^ 
frame, i.e., Eqs. (|^) and (^, even though the nonvalence contribution in Fig. 3(c) is absent 
here. 

After some manipulation, we obtain the EM form factor of a qQ bound state in the 
g"*" = frame as follows 

jr(g2^g+ = 0) = ^ / deil-a-hsmO) 

^q^{q^ - 4M2) Arcsin{^f ) 



X tanh ^ 



q^ a + h sin 9 



(23) 



g2 - 4M2 ih cos 6^ 

where a = (x+ + x„)/2, b = (a;+ — x_)/2, and g^ = — g^. While the representations in 
Eqs. (|T^) and (^) look apparently different, the two formulas, Eqs. (|T^ and (PBD, turn 



out to be actually identical. As we will show explicitly in the next section of numerical 
calculations, all three results of Eqs. ([10|) , (|T9D and (|23|) indeed coincide exactly in the 
entire g^ range. 

IV. NUMERICAL RESULTS 

For our numerical analysis of n, K, and D meson form factors, we use the physical meson 
masses together with the following constituent quark and antiquark masses: m„ = m^ = 0.25 
GeV, rris = 0.48 GeV, and nic = 1.8 GeV P,[7|j21|l. Since our numerical results of the EM 
form factors obtained from Eqs. {^^, ([19|) and ( p3D turn out to be exactly same with each 
other for the entire g^ region, only a single line is depicted in Figs. 4, 5 and 6 for the form 
factor calculations of tt, i^, and D mesons, respectively. 

In should be noted from our constituent masses that M^ < m'^ + m% for vr and K and 
M^ > rrig + m% for D meson cases. As discussed in Ref. [|18[ for the analysis of the one- 
particle matrix element of a scalar current, the sigularity for M^ > m"^ + m% case starts 
at 
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Cin = Z;^K(Q) - (^ - "^QW) ][(^ + ^Q(,)) - S(Q)] (24) 



for 'y*q(j{'y*QQ) vertex, while the singularity for M^ < m"^ + ?7i?j case starts on the positive 
g^-axis at the threshold point g^j^ = 4m^,g^ for 'y*qq{'y*QQ) vertex. Our numerical results 
exhibit all of these threshold behaviors coming from the normal (vr, K) and anomalous (D) 
cases. As a consistency check, we also compare our numerical results of the form factor 
F{q^) = Re F{q^) + z Im F{q^) with the dispersion relations given by 



Reftf) = iprlHL^,,", (25) 

-K J-oD q'^ — q^ 

\^F{q^) = -^-Pr^^^f^dq'\ (26) 



where P indicates the Cauchy principal value. 

In Fig. 4(a), we show the EM form factor of the pion for —2 GeV^ < g^ < 3 GeV^. The 
imaginary part (the dotted line) of the form factor starts at g^j^ = 4mL^) = 0.25 GeV, which 
is consistent with the condition for M^ < m? + tuq case. It is interesting to note that the 
square of the total form factor |F^(g^)p (thick solid line) produces a p meson-type peak near 
g2 r^ ]^^2 However, we do not claim that this model indeed reproduces all the features of the 
vector meson dominance (VMD) phenomena because the more realistic phenomenological 
models may have to incorporate the more complex mechanism including the final state 
interaction. More investigation along this line is under consideration. Nevertheless, this 
simple model is capable of generating the peak and the position of peak is quite consistent 
with the VMD. 

In Fig. 4(b), we show the timelike form factor of the pion for the entire g^ > re- 
gion and compare the imaginary part of our direct calculations (dotted line) obtained from 
Eqs. (p!OD, (p^D, and (^) with the result (data of black dots) obtained from the dispersion 
relations given by Eq. (|26D. Our direct calculation is in an excellent agreement with the 
solution of the dispersion relations. Our results for the real part are also confirmed to be in 
complete agreement with the dispersion relations. For high g^ region, the imaginary part of 
the form factor is dominant over the real part (thin solid line). 
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In Fig. 5(a), we show the kaon form factor for —2 GeV^ < Q'^ < 5 GeV^. The kaon also 
has the normal singularity. However, it has two thresholds for the imaginary parts; one is 
?min — 4"^M ^^^ ^h^ other is g^j^ = 4m^. These lead to the humped shape (dotted line) of 
the imaginary part shown in Fig. 5(a). While we have in principle two vector- meson-type 
peaks (i.e. p and 0), one can see in Fig. 5(a) only meson-type peak for the timelike kaon 
EM form factor above the physical threshold at g^j^ = 4M]^. We also show in Fig. 5(b) 
the imaginary part from our direct calculation is in an excellent agreement with the result 
(data of black dots) from the dispersion relations for the entire timelike q"^ region. Again, 
the imaginary part is predominant for high g^ region. 

In Fig. 6, we show the D meson form factor for —10 GeV^ < g^ < 30 GeV^. Unlike 
the normal threshold of tt and K form factor calculations, the D meson form factor shows 
anomalous thresholds according to Eq. (^4]) , i.e., g^j^ ~ 0.24 GeV^ (compared to 4m^ = 0.25 
GeV^ for normal case) and g^jj^ ~ 12.4 GeV^ (compared to 4m^ = 12.96 GeV^ for normal 
case) for the 7* — c? and 7* — c vertices, respectively. Similar to the kaon case in Fig. 5, we 
also have two unphysical peaks, i.e., p and J /ip{lS) meson type peaks due to J and c quarks, 
respectively. However, the timelike form factor of D meson has no pole structure for the the 
physical g^ > 4M|, region. In all of these figures (Figs. 4-6), the numerical result of Eq. ( |23D 
obtained from g^ = frame without encountering the nonvalence diagram coincides exactly 
with the numerical results of Eqs. (jT^) and (|TU|) obtained from g"*" 7^ frame. 

V. CONCLUSION AND DISCUSSION 

While it is not a simple task to explore the timelike region with the framework at hand, 
we attempted in this work to provide a clear example demonstrating that one can compute 
the timelike form factor without encountering the nonvalence contributions. We investigated 
both spacelike and timelike regions of the meson EM form factor using an exactly solvable 
model of (3 + 1) dimensional scalar field theory interacting with gauge fields (see Eq. (|1])). 
We calculated the form factor in spacelike region using the Drell-Yan-West (g"*" = 0) frame 
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and showed that its analytic continuation to the timehke region reproduces exactly the 
direct result of timelike form factor obtained in the longitudinal momentum (g+ 7^ and 
q± = 0) frame. The analytic continuation of the result in the g+ = frame to the timelike 
region automatically generates the effect of the nonvalence contributions to the timelike form 
factor. Another interesting result in our model calculations is that the peaks analogous to 
the VMD were generated and the position of peaks were indeed quite consistent with the 
VMD. Even though phenomenologically more elaborate analyses including the final state 
interaction may be necessary to reproduce the more realistic feature of VMD, our results 
seem pretty encouraging for further investigations. Using the dispersion relations, we have 
also confirmed that our numerical results of imaginary parts start at q'^^^ = 4mg(Q) and 
the normal thresholds appear for vr and K (M^ < m'^ + 171%) systems while the anomalous 
threshold exists for D (M^ > m'^ + 171%) system. However, we note that our form factor 
calculations are effectively various ways of evaluating the Feynman perturbation theory 
triangle diagrams and thus the analytic properties presented in sections II, III, and IV are 
restricted to one-loop diagrams. 

Even though the model (Eq. (Q)) may not be as realistic as one may want, we think 
that it plays the role of some guidance for the more realistic phenomenological models. For 
an application of our analytic continuation method to the more realistic model, we have 
considered an immediate extension of bosonic contact interaction to the fermionic covariant 
model with the monopole smearing vertices used by de Melo et al. [0,0 and investigated 
singular free processes in timelike region such as 0^(= J^) -^ 0^ and 0^ -^ 1^ semileptonic 
decays. In this application, we found that the covariant solution of the weak form factors 
/+ for 0^ -^ 0^ and V and A2 for 0^ -^ 1^ processes are exactly equal to our analytic 
solutions obtained from the current j+ in g+ = frame where only valence diagram is 
needed. On the other hand, the form factors /_ for 0^ -^ 0~ and Aq and Ai for 0^ ^ 1~ 
processes turn out to be affected by the zero-mode contributions and one has to include the 
zero-mode effects to make our analytic solutions from g+ = frame equal to the covariant 
solutions (see also |]10|)- This covariant model calculation further verifies the validity of 
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our analytic continuation method. We also applied our analytic continuation method to 
the same semileptonic decay processes of various pseudoscalar and vector mesons using our 
LFQM []7|, |2T] , |22| and found a good agreement with the available experimental data as well 
as the lattice calculations. 

However, an additional step of identifying the zero-mode contributions is necessary in 
handling the zero-mode issue associated with the currents j± for 0~ -^ 0~ and j~^ coupled 
with longitudinal decay modes for 0~ -^ 1' exclusive semileptonic decay processes. Also, 
the singular structures around the pole-regions in the timelike elastic form factors needs to 
be further investigated in the more phenomenological model. Nevertheless, it is encouraging 
that in this first attempt of calculating the timelike elastic meson form factor, the analytic 
continuation method seems to work as nicely as one may have anticipated. 
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FIG. 4. a. The electromagnetic form fator of the pion in (3 + 1) dimensional scalar field theory 
for —2 < g^ < 3 GeV^. The total, real, and imaginary parts of \F.,^{q^)\'^ are represented by thick 
solid, solid, and dotted lines, respectively. 
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FIG. 4.b. The electromagnetic form fator of the pion in (3 + 1) dimensional scalar field theory 
for the entire timelike region compared to the dispersion relations (data of black dots) given by 
Eq. (25). The same line code as in Fig. 4(a) is used. 
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FIG. 5. a. The electromagnetic form fator of the kaon in (3 + 1) dimensional scalar field theory 
for —2 < g^ < 5 GeV^. The same line code as in Fig. 4(a) is used. 
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FIG. 5.b. The electromagnetic form fator of the kaon in (3 + 1) dimensional scalar field theory 
for the entire timelike region compared to the dispersion relations (data of black dots) given by 
Eq. (25). The same line code as in Fig. 4(a) is used. 
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FIG. 6. The electromagnetic form fator of the D meson in (3 + 1) dimensional scalar field theory 
for — 10 < (7^ < 30 GeV^. The same line code as in Fig. 4(a) is used. 
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